Viruses associated with measles-like illnesses in Uganda.

OBJECTIVES: In this study, we investigated the causes of measles-like illnesses (MLI) in the Uganda national surveillance program in order to inform diagnostic assay selection and vaccination strategies. METHODS: We used metagenomic next-generation sequencing (M-NGS) on the Illumina platform to identify viruses associated with MLI (defined as fever and rash in the presence of either cough, coryza or conjunctivitis) in patient samples that had tested IgM negative for measles between 2010 and 2019. RESULTS: Viral genomes were identified in 87/271 (32%) of samples, of which 44/271 (16%) contained 12 known viral pathogens. Expected viruses included rubella, human parvovirus B19, Epstein Barr virus, human herpesvirus 6B, human cytomegalovirus, varicella zoster virus and measles virus (detected within the seronegative window-period of infection) and the blood-borne hepatitis B virus. We also detected Saffold virus, human parvovirus type 4, the human adenovirus C2 and vaccine-associated poliovirus type 1. CONCLUSIONS: The study highlights the presence of undiagnosed viruses causing MLI in Uganda, including vaccine-preventable illnesses. NGS can be used to monitor common viral infections at a population level, especially in regions where such infections are prevalent, including low and middle income countries to guide vaccination policy and optimize diagnostic assays.

Adeno-associated virus 2 infection in children with non-A-E hepatitis.

An outbreak of acute hepatitis of unknown aetiology in children was reported in Scotland1 in April 2022 and has now been identified in 35 countries2. Several recent studies have suggested an association with human adenovirus with this outbreak, a virus not commonly associated with hepatitis. Here we report a detailed case-control investigation and find an association between adeno-associated virus 2 (AAV2) infection and host genetics in disease susceptibility. Using next-generation sequencing, PCR with reverse transcription, serology and in situ hybridization, we detected recent infection with AAV2 in plasma and liver samples in 26 out of 32 (81%) cases of hepatitis compared with 5 out of 74 (7%) of samples from unaffected individuals. Furthermore, AAV2 was detected within ballooned hepatocytes alongside a prominent T cell infiltrate in liver biopsy samples. In keeping with a CD4+ T-cell-mediated immune pathology, the human leukocyte antigen (HLA) class II HLA-DRB1*04:01 allele was identified in 25 out of 27 cases (93%) compared with a background frequency of 10 out of 64 (16%; P = 5.49 × 10-12). In summary, we report an outbreak of acute paediatric hepatitis associated with AAV2 infection (most likely acquired as a co-infection with human adenovirus that is usually required as a 'helper virus' to support AAV2 replication) and disease susceptibility related to HLA class II status.

Genomic and Phenotypic Comparisons Reveal Distinct Variants of Wolbachia Strain wAlbB.

The intracellular bacterium Wolbachia inhibits virus replication and is being harnessed around the world to fight mosquito-borne diseases through releases of mosquitoes carrying the symbiont. Wolbachia strains vary in their ability to invade mosquito populations and suppress viruses in part due to differences in their density within the insect and associated fitness costs. Using whole-genome sequencing, we demonstrate the existence of two variants in wAlbB, a Wolbachia strain being released in natural populations of Aedes aegypti mosquitoes. The two variants display striking differences in genome architecture and gene content. Differences in the presence/absence of 52 genes between variants include genes located in prophage regions and others potentially involved in controlling the symbiont's density. Importantly, we show that these genetic differences correlate with variation in wAlbB density and its tolerance to heat stress, suggesting that different wAlbB variants may be better suited for field deployment depending on local environmental conditions. Finally, we found that the wAlbB genome remained stable following its introduction in a Malaysian mosquito population. Our results highlight the need for further genomic and phenotypic characterization of Wolbachia strains in order to inform ongoing Wolbachia-based programs and improve the selection of optimal strains in future field interventions. IMPORTANCE Dengue is a viral disease transmitted by Aedes mosquitoes that threatens around half of the world population. Recent advances in dengue control involve the introduction of Wolbachia bacterial symbionts with antiviral properties into mosquito populations, which can lead to dramatic decreases in the incidence of the disease. In light of these promising results, there is a crucial need to better understand the factors affecting the success of such strategies, in particular the choice of Wolbachia strain for field releases and the potential for evolutionary changes. Here, we characterized two variants of a Wolbachia strain used for dengue control that differ at the genomic level and in their ability to replicate within the mosquito. We also found no evidence for the evolution of the symbiont within the 2 years following its deployment in Malaysia. Our results have implications for current and future Wolbachia-based health interventions.

The apparent interferon resistance of transmitted HIV-1 is possibly a consequence of enhanced replicative fitness.

HIV-1 transmission via sexual exposure is an inefficient process. When transmission does occur, newly infected individuals are colonized by the descendants of either a single virion or a very small number of establishing virions. These transmitted founder (TF) viruses are more interferon (IFN)-resistant than chronic control (CC) viruses present 6 months after transmission. To identify the specific molecular defences that make CC viruses more susceptible to the IFN-induced 'antiviral state', we established a single pair of fluorescent TF and CC viruses and used arrayed interferon-stimulated gene (ISG) expression screening to identify candidate antiviral effectors. However, we observed a relatively uniform ISG resistance of transmitted HIV-1, and this directed us to investigate possible underlying mechanisms. Simple simulations, where we varied a single parameter, illustrated that reduced growth rate could possibly underly apparent interferon sensitivity. To examine this possibility, we closely monitored in vitro propagation of a model TF/CC pair (closely matched in replicative fitness) over a targeted range of IFN concentrations. Fitting standard four-parameter logistic growth models, in which experimental variables were regressed against growth rate and carrying capacity, to our in vitro growth curves, further highlighted that small differences in replicative growth rates could recapitulate our in vitro observations. We reasoned that if growth rate underlies apparent interferon resistance, transmitted HIV-1 would be similarly resistant to any growth rate inhibitor. Accordingly, we show that two transmitted founder HIV-1 viruses are relatively resistant to antiretroviral drugs, while their matched chronic control viruses were more sensitive. We propose that, when present, the apparent IFN resistance of transmitted HIV-1 could possibly be explained by enhanced replicative fitness, as opposed to specific resistance to individual IFN-induced defences. However, further work is required to establish how generalisable this mechanism of relative IFN resistance might be.

Genome sequencing and comparative analysis of Wolbachia strain wAlbA reveals Wolbachia-associated plasmids are common.

Wolbachia are widespread maternally-transmitted bacteria of arthropods that often spread by manipulating their host's reproduction through cytoplasmic incompatibility (CI). Their invasive potential is currently being harnessed in field trials aiming to control mosquito-borne diseases. Wolbachia genomes commonly harbour prophage regions encoding the cif genes which confer their ability to induce CI. Recently, a plasmid-like element was discovered in wPip, a Wolbachia strain infecting Culex mosquitoes; however, it is unclear how common such extra-chromosomal elements are in Wolbachia. Here we sequenced the complete genome of wAlbA, a strain of the symbiont found in Aedes albopictus, after eliminating the co-infecting and higher density wAlbB strain that previously made sequencing of wAlbA challenging. We show that wAlbA is associated with two new plasmids and identified additional Wolbachia plasmids and related chromosomal islands in over 20% of publicly available Wolbachia genome datasets. These plasmids encode a variety of accessory genes, including several phage-like DNA packaging genes as well as genes potentially contributing to host-symbiont interactions. In particular, we recovered divergent homologues of the cif genes in both Wolbachia- and Rickettsia-associated plasmids. Our results indicate that plasmids are common in Wolbachia and raise fundamental questions around their role in symbiosis. In addition, our comparative analysis provides useful information for the future development of genetic tools to manipulate and study Wolbachia symbionts.

Characterisation of a Hepatitis C Virus Subtype 2a Cluster in Scottish PWID with a Suboptimal Response to Glecaprevir/Pibrentasvir Treatment.

Direct-acting antivirals (DAAs) have revolutionised the treatment of Hepatitis C virus (HCV), allowing the World Health Organisation (WHO) to set a target of eliminating HCV by 2030. In this study we aimed to investigate glecaprevir and pibrentasvir (GP) treatment outcomes in a cohort of patients with genotype 2a infection. METHODS: Clinical data and plasma samples were collected in NHS Greater Glasgow & Clyde. Next generation whole genome sequencing and replicon assays were carried out at the MRC-University of Glasgow Centre for Virus Research. RESULTS: 132 cases infected with genotype 2a HCV were identified. The SVR rate for this group was 91% (112/123) following treatment with GP. An NS5A polymorphism, L31M, was detected in all cases of g2a infection, and L31M+R353K in individuals that failed treatment. The results showed that R353K was present in 90% of individuals in the Glasgow genotype 2a phylogenetic cluster but in less than 5% of all HCV subtype 2a published sequences. In vitro efficacy of pibrentasvir against sub-genomic replicon constructs containing these mutations showed a 2-fold increase in IC50 compared to wildtype. CONCLUSION: This study describes a cluster of HCV genotype 2a infection associated with a lower-than-expected SVR rate following GP treatment in association with the NS5A mutations L31M+R353K.

Genomic epidemiology of SARS-CoV-2 in a university outbreak setting and implications for public health planning.

Whole genome sequencing of SARS-CoV-2 has occurred at an unprecedented scale, and can be exploited for characterising outbreak risks at the fine-scale needed to inform control strategies. One setting at continued risk of COVID-19 outbreaks are higher education institutions, associated with student movements at the start of term, close living conditions within residential halls, and high social contact rates. Here we analysed SARS-CoV-2 whole genome sequences in combination with epidemiological data to investigate a large cluster of student cases associated with University of Glasgow accommodation in autumn 2020, Scotland. We identified 519 student cases of SARS-CoV-2 infection associated with this large cluster through contact tracing data, with 30% sequencing coverage for further analysis. We estimated at least 11 independent introductions of SARS-CoV-2 into the student population, with four comprising the majority of detected cases and consistent with separate outbreaks. These four outbreaks were curtailed within a week following implementation of control measures. The impact of student infections on the local community was short-term despite an underlying increase in community infections. Our study highlights the need for context-specific information in the formation of public health policy for higher educational settings.

SARS-CoV-2 Omicron is an immune escape variant with an altered cell entry pathway.

Vaccines based on the spike protein of SARS-CoV-2 are a cornerstone of the public health response to COVID-19. The emergence of hypermutated, increasingly transmissible variants of concern (VOCs) threaten this strategy. Omicron (B.1.1.529), the fifth VOC to be described, harbours multiple amino acid mutations in spike, half of which lie within the receptor-binding domain. Here we demonstrate substantial evasion of neutralization by Omicron BA.1 and BA.2 variants in vitro using sera from individuals vaccinated with ChAdOx1, BNT162b2 and mRNA-1273. These data were mirrored by a substantial reduction in real-world vaccine effectiveness that was partially restored by booster vaccination. The Omicron variants BA.1 and BA.2 did not induce cell syncytia in vitro and favoured a TMPRSS2-independent endosomal entry pathway, these phenotypes mapping to distinct regions of the spike protein. Impaired cell fusion was determined by the receptor-binding domain, while endosomal entry mapped to the S2 domain. Such marked changes in antigenicity and replicative biology may underlie the rapid global spread and altered pathogenicity of the Omicron variant.

The evolutionary dynamics and epidemiological history of hepatitis C virus genotype 6, including unique strains from the Li community of Hainan Island, China.

Hepatitis C virus (HCV) is a highly diverse pathogen that frequently establishes a chronic long-term infection, but the origins and drivers of HCV diversity in the human population remain unclear. Previously unidentified strains of HCV genotype 6 (gt6) were recently discovered in chronically infected individuals of the Li ethnic group living in Baisha County, Hainan Island, China. The Li community, who were early settlers on Hainan Island, has a distinct host genetic background and cultural identity compared to other ethnic groups on the island and mainland China. In this report, we generated 33 whole virus genome sequences to conduct a comprehensive molecular epidemiological analysis of these novel gt6 strains in the context of gt6 isolates present in Southeast Asia. With the exception of one gt6a isolate, the Li gt6 sequences formed three novel clades from two lineages which constituted 3 newly assigned gt6 subtypes and 30 unassigned strains. Using Bayesian inference methods, we dated the most recent common ancestor for all available gt6 whole virus genome sequences to approximately 2767 bce (95 per cent highest posterior density (HPD) intervals, 3670-1397 bce), which is far earlier than previous estimates. The substitution rate was 1.20 × 10-4 substitutions/site/year (s/s/y), and this rate varied across the genome regions, from 1.02 × 10-5 s/s/y in the 5'untranslated region (UTR) region to 3.07 × 10-4 s/s/y in E2. Thus, our study on an isolated ethnic minority group within a small geographical area of Hainan Island has substantially increased the known diversity of HCV gt6, already acknowledged as the most diverse HCV genotype. The extant HCV gt6 sequences from this study were probably transmitted to the Li through at least three independent events dating perhaps from around 4,000 years ago. This analysis describes deeper insight into basic aspects of HCV gt6 molecular evolution including the extensive diversity of gt6 sequences in the isolated Li ethnic group.

A clinical and molecular epidemiological survey of hepatitis C in Blantyre, Malawi, suggests a historic mechanism of transmission.

Hepatitis C virus (HCV) is a leading cause of liver disease worldwide. There are no previous representative community HCV prevalence studies from Southern Africa, and limited genotypic data. Epidemiological data are required to inform an effective public health response. We conducted a household census-based random sampling serological survey, and a prospective hospital-based study of patients with cirrhosis and hepatocellular carcinoma (HCC) in Blantyre, Malawi. We tested participants with an HCV antigen/antibody ELISA (Monolisa, Bio-Rad), confirmed with PCR (GeneXpert, Cepheid) and used line immunoassay (Inno-LIA, Fujiribio) for RNA-negative participants. We did target-enrichment whole-genome HCV sequencing (NextSeq, Illumina). Among 96,386 censused individuals, we randomly selected 1661 people aged ≥16 years. Population-standardized HCV RNA prevalence was 0.2% (95% CI 0.1-0.5). Among 236 patients with cirrhosis and HCC, HCV RNA prevalence was 1.9% and 5.0%, respectively. Mapping showed that HCV RNA+ patients were from peri-urban areas surrounding Blantyre. Community and hospital HCV RNA+ participants were older than comparator HCV RNA-negative populations (median 53 vs 30 years for community, p = 0.01 and 68 vs 40 years for cirrhosis/HCC, p < 0.001). Endemic HCV genotypes (n = 10) were 4v (50%), 4r (30%) and 4w (10%). In this first census-based community serological study in Southern Africa, HCV was uncommon in the general population, was centred on peri-urban regions and was attributable for <5% of liver disease. HCV infection was observed only among older people, suggesting a historic mechanism of transmission. Genotype 4r, which has been associated with treatment failure with ledipasvir and daclatasvir, is endemic.

Real-World Outcomes of Direct-Acting Antiviral Treatment and Retreatment in United Kingdom-Based Patients Infected With Hepatitis C Virus Genotypes/Subtypes Endemic in Africa.

BACKGROUND: Chronic hepatitis C virus (HCV) infection affects 71 million individuals, mostly residing in low- and middle-income countries (LMICs). Direct-acting antivirals (DAAs) give high rates of sustained virological response (SVR) in high-income countries where a restricted range of HCV genotypes/subtypes circulate. METHODS: We studied United Kingdom-resident patients born in Africa to examine DAA effectiveness in LMICs where there is far greater breadth of HCV genotypes/subtypes. Viral genome sequences were determined from 233 patients. RESULTS: Full-length viral genomic sequences for 26 known subtypes and 5 previously unidentified isolates covering 5 HCV genotypes were determined. From 149 patients who received DAA treatment/retreatment, the overall SVR was 93%. Treatment failure was associated primarily with 2 subtypes, gt1l and gt4r, using sofosbuvir/ledipasvir. These subtypes contain natural resistance-associated variants that likely contribute to poor efficacy with this drug combination. Treatment failure was also significantly associated with hepatocellular carcinoma. CONCLUSIONS: DAA combinations give high SVR rates despite the high HCV diversity across the African continent except for subtypes gt1l and gt4r, which respond poorly to sofosbuvir/ledipasvir. These subtypes are widely distributed across Western, Central, and Eastern Africa. Thus, in circumstances where accurate genotyping is absent, ledipasvir and its generic compounds should not be considered as a recommended treatment option.

Long-term adaptation following influenza A virus host shifts results in increased within-host viral fitness due to higher replication rates, broader dissemination within the respiratory epithelium and reduced tissue damage.

The mechanisms and consequences of genome evolution on viral fitness following host shifts are poorly understood. In addition, viral fitness -the ability of an organism to reproduce and survive- is multifactorial and thus difficult to quantify. Influenza A viruses (IAVs) circulate broadly among wild birds and have jumped into and become endemic in multiple mammalian hosts, including humans, pigs, dogs, seals, and horses. H3N8 equine influenza virus (EIV) is an endemic virus of horses that originated in birds and has been circulating uninterruptedly in equine populations since the early 1960s. Here, we used EIV to quantify changes in infection phenotype associated to viral fitness due to genome-wide changes acquired during long-term adaptation. We performed experimental infections of two mammalian cell lines and equine tracheal explants using the earliest H3N8 EIV isolated (A/equine/Uruguay/63 [EIV/63]), and A/equine/Ohio/2003 (EIV/2003), a monophyletic descendant of EIV/63 isolated 40 years after the emergence of H3N8 EIV. We show that EIV/2003 exhibits increased resistance to interferon, enhanced viral replication, and a more efficient cell-to-cell spread in cells and tissues. Transcriptomics analyses revealed virus-specific responses to each virus, mainly affecting host immunity and inflammation. Image analyses of infected equine respiratory explants showed that despite replicating at higher levels and spreading over larger areas of the respiratory epithelium, EIV/2003 induced milder lesions compared to EIV/63, suggesting that adaptation led to reduced tissue pathogenicity. Our results reveal previously unknown links between virus genotype and the host response to infection, providing new insights on the relationship between virus evolution and fitness.

Genetic epidemiology of SARS-CoV-2 transmission in renal dialysis units - A high risk community-hospital interface.

OBJECTIVES: Patients requiring haemodialysis are at increased risk of serious illness with SARS-CoV-2 infection. To improve the understanding of transmission risks in six Scottish renal dialysis units, we utilised the rapid whole-genome sequencing data generated by the COG-UK consortium. METHODS: We combined geographical, temporal and genomic sequence data from the community and hospital to estimate the probability of infection originating from within the dialysis unit, the hospital or the community using Bayesian statistical modelling and compared these results to the details of epidemiological investigations. RESULTS: Of 671 patients, 60 (8.9%) became infected with SARS-CoV-2, of whom 16 (27%) died. Within-unit and community transmission were both evident and an instance of transmission from the wider hospital setting was also demonstrated. CONCLUSIONS: Near-real-time SARS-CoV-2 sequencing data can facilitate tailored infection prevention and control measures, which can be targeted at reducing risk in these settings.

Subgenomic RNA identification in SARS-CoV-2 genomic sequencing data.

We have developed periscope, a tool for the detection and quantification of subgenomic RNA (sgRNA) in SARS-CoV-2 genomic sequence data. The translation of the SARS-CoV-2 RNA genome for most open reading frames (ORFs) occurs via RNA intermediates termed "subgenomic RNAs." sgRNAs are produced through discontinuous transcription, which relies on homology between transcription regulatory sequences (TRS-B) upstream of the ORF start codons and that of the TRS-L, which is located in the 5' UTR. TRS-L is immediately preceded by a leader sequence. This leader sequence is therefore found at the 5' end of all sgRNA. We applied periscope to 1155 SARS-CoV-2 genomes from Sheffield, United Kingdom, and validated our findings using orthogonal data sets and in vitro cell systems. By using a simple local alignment to detect reads that contain the leader sequence, we were able to identify and quantify reads arising from canonical and noncanonical sgRNA. We were able to detect all canonical sgRNAs at the expected abundances, with the exception of ORF10. A number of recurrent noncanonical sgRNAs are detected. We show that the results are reproducible using technical replicates and determine the optimum number of reads for sgRNA analysis. In VeroE6 ACE2+/- cell lines, periscope can detect the changes in the kinetics of sgRNA in orthogonal sequencing data sets. Finally, variants found in genomic RNA are transmitted to sgRNAs with high fidelity in most cases. This tool can be applied to all sequenced COVID-19 samples worldwide to provide comprehensive analysis of SARS-CoV-2 sgRNA.

Circulating SARS-CoV-2 spike N439K variants maintain fitness while evading antibody-mediated immunity.

SARS-CoV-2 can mutate and evade immunity, with consequences for efficacy of emerging vaccines and antibody therapeutics. Here, we demonstrate that the immunodominant SARS-CoV-2 spike (S) receptor binding motif (RBM) is a highly variable region of S and provide epidemiological, clinical, and molecular characterization of a prevalent, sentinel RBM mutation, N439K. We demonstrate N439K S protein has enhanced binding affinity to the hACE2 receptor, and N439K viruses have similar in vitro replication fitness and cause infections with similar clinical outcomes as compared to wild type. We show the N439K mutation confers resistance against several neutralizing monoclonal antibodies, including one authorized for emergency use by the US Food and Drug Administration (FDA), and reduces the activity of some polyclonal sera from persons recovered from infection. Immune evasion mutations that maintain virulence and fitness such as N439K can emerge within SARS-CoV-2 S, highlighting the need for ongoing molecular surveillance to guide development and usage of vaccines and therapeutics.

Genomic epidemiology reveals multiple introductions of SARS-CoV-2 from mainland Europe into Scotland.

Coronavirus disease 2019 (COVID-19) was first diagnosed in Scotland on 1 March 2020. During the first month of the outbreak, 2,641 cases of COVID-19 led to 1,832 hospital admissions, 207 intensive care admissions and 126 deaths. We aimed to identify the source and number of introductions of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) into Scotland using a combined phylogenetic and epidemiological approach. Sequencing of 1,314 SARS-CoV-2 viral genomes from available patient samples enabled us to estimate that SARS-CoV-2 was introduced to Scotland on at least 283 occasions during February and March 2020. Epidemiological analysis confirmed that early introductions of SARS-CoV-2 originated from mainland Europe (the majority from Italy and Spain). We identified subsequent early outbreaks in the community, within healthcare facilities and at an international conference. Community transmission occurred after 2 March, 3 weeks before control measures were introduced. Earlier travel restrictions or quarantine measures, both locally and internationally, would have reduced the number of COVID-19 cases in Scotland. The risk of multiple reintroduction events in future waves of infection remains high in the absence of population immunity.

Metaviromics Reveals Unknown Viral Diversity in the Biting Midge Culicoides impunctatus.

Biting midges (Culicoides species) are vectors of arboviruses and were responsible for the emergence and spread of Schmallenberg virus (SBV) in Europe in 2011 and are likely to be involved in the emergence of other arboviruses in Europe. Improved surveillance and better understanding of risks require a better understanding of the circulating viral diversity in these biting insects. In this study, we expand the sequence space of RNA viruses by identifying a number of novel RNA viruses from Culicoides impunctatus (biting midge) using a meta-transcriptomic approach. A novel metaviromic pipeline called MetaViC was developed specifically to identify novel virus sequence signatures from high throughput sequencing (HTS) datasets in the absence of a known host genome. MetaViC is a protein centric pipeline that looks for specific protein signatures in the reads and contigs generated as part of the pipeline. Several novel viruses, including an alphanodavirus with both segments, a novel relative of the Hubei sobemo-like virus 49, two rhabdo-like viruses and a chuvirus, were identified in the Scottish midge samples. The newly identified viruses were found to be phylogenetically distinct to those previous known. These findings expand our current knowledge of viral diversity in arthropods and especially in these understudied disease vectors.